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Abstract. We describe a general method of arithmetic coding of geodesies 
on the modular surface based on the study of one-dimensional Gauss-like maps 
associated to a two parameter family of continued fractions introduced in 1161. 
The finite rectangular structure of the attractors of the natural extension maps 
and the corresponding "reduction theory" play an essential role. In special 
cases, when an (a, 6)-expansion admits a so-called "dual" , the coding sequences 
are obtained by juxtaposition of the boundary expansions of the fixed points, 
and the set of coding sequences is a countable sofic shift. We also prove that 
the natural extension maps are Bernoulli shifts and compute the density of the 
absolutely continuous invariant measure and the measure-theoretic entropy of 
the one-dimensional map. 



1. Introduction and background 



In [16], the authors studied a new two- parameter family of continued fraction 
transformations. These transformations can be defined using the standard gener- 
ators T{x) = X + 1, S{x) = —1/x of the modular group 5^(2, Z) and considering 
: t ^ t given by 



X -\- 1 if X < a 
(1-1) fa,b{^)=\-- ifa<x<b 

X 

X — 1 if X > b . 

Under the assumption that the parameters (a, b) belong to the set 

V = {(a, b)\a <0 <b, b- a>l, -ab <!} , 

one can introduce corresponding continued fraction algorithms by using the first 
return map of fa^b to the interval [a, b). Equivalently, these so called (a, b) -continued 
fractions can be defined using a generalized integral part function: 



(1.2) 



— aj if X < a 
if a <x <b 

Jx — b] if X >b , 

where [x\ denotes the integer part of x and [x] = [x\ + 1. 
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A starting point of the theory is the fohowing result [16 , Theorem 2.1]: if (a, b) G 
7^, then any irrational number x can be expressed uniquely as an infinite continued 
fraction of the form 

X = no = [no.ni, • • •]a,6, (^fc 7^ for > 1), 

ni 

n2 

where no = L^la,6, ^1 = -^3^ and n/e+i = L^fe+ila,6, ^fe+i = - xk-uk ' 
sequence of partial fractions rk = [no, ni, . . . , nk']a,b converges to x. 

It is possible to construct (a, 6) -continued fraction expansions for rational num- 
bers, too. However, such expansions will terminate after finitely many steps if 
b ^ 0. If 6 = 0, the expansions of rational numbers will end with a tail of 2's, since 
0= [1,2,2,. ..la,o. 

The above family of continued fraction transformations contains three classical 
examples: the case a = — 1, 6 = described in [22l [TJ gives the "minus" (back- 
ward) continued fractions, the case a = —1/2, b = 1/2 gives the "closest-integer" 
continued fractions considered first by Hurwitz in [9 , and the case a = — 1, 6 = 1 
was presented in [W, 14 in connection with a method of coding symbolically the 
geodesic fiow on the modular surface following Artin's pioneering work [6 and cor- 
responds to the regular "plus" continued fractions with alternating signs of the 
digits. 

The main object of study in [16 is a two-dimensional realization of the natural 
extension map of /a, 5, Fa, 6 : \ A ^ \ A, A = {(x, y) G M^|x = y}, defined by 

{(x + 1,^ + 1) if ^<a 
(x- 1,^-1) ify>b. 
Here is the main result of that paper: 

Theorem 1.1 ([16 ). There exists an explicit one- dimensional Lebesgue measure 
zero, uncountable set £ that lies on the diagonal boundary b = a-\-l ofV such that: 

(1) for all (a, b) eV\£ the map Fa^h has an attractor Da^b = n^_o^a,6(^^ \^) 
on which Fa^h ^-^ essentially bijective. 

(2) The set Da^b consists of two (or one, in degenerate cases) connected com- 
ponents each having finite rectangular structure^ i.e. bounded by non- 
decreasing step -functions with a finite number of steps. 

(3) Almost every point {x^y) of the plane (x 7^ y) is mapped to Da^b after 
finitely many iterations of Fa^b- 

An essential role in the argument is played by the forward orbits associated to a 
and b: to a, the upper orbit Ou{o) (i.e. the orbit of So) and the lower orbit Oi{a) 
(i.e. the orbit of Ta), and to 6, the upper orbit Ou{b) (i.e. the orbit of T~^b) and the 
lower orbit Oi{b) (i.e. the orbit of Sb). It was proved in [16] that if (a, b) e V\S^ 
then fa^b satisfies the finiteness condition, i.e. for both a and 6, their upper and 
lower orbits are either eventually periodic, or they satisfy the cycle property, i.e. 
they meet forming a cycle; more precisely, there exist /ci,mi, /c2,m2 > s.t. 
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Figure 1. Attracting domain Da^b for a = - |, 6 = | 



where Ca and C5 are the ends of the cycles. If the products of transformations 
over the upper and lower sides of the cycle are equal, the cycle property is strong^ 
otherwise, it is weak. In both cases the set Ca,b of the corresponding values is finite; 
ends of the cycles belong to the set £^,6 if and only if they are equal to 0, i.e. if the 
cycle is weak. The structure of the attractor Da^b is explicitly "computed" from 
the finite set Ca^b- 

The paper is organized as follows. In Section 2 we give some background infor- 
mation about geodesic flows and their representations as special flows over symbolic 
dynamical systems, and deflne the coding map. In Section 3 we describe the reduc- 
tion procedure for coding geodesies via (a, 6)-continued fractions based on the study 
of the attractor of the associated natural extension map, deflne the corresponding 
cross-section set, and introduce the notion of reduced geodesic. In Section 4 we prove 
that the flrst return map to the cross-section corresponds to a shift of the coding 
sequence (Theorem |4.1| ) and, as a consequence, show that (a, 6)-continued fractions 
satisfy the Tail Property, i.e. two 51/(2, Z)-equivalent real numbers have the same 
tails in their (a, 6)-continued fraction expansions. In Section 5 we introduce a notion 
of a dual code and show that if an (a, 6)-expansion has a dual (a', 6')-expansion, then 
the coding sequence of a reduced geodesic is obtained by juxtaposition of the (a, b)- 
expansion of its attracting endpoint w and the (a^, 6')-expansion of where u is 
its repelling endpoint. We also prove that if the (a, 6)-expansion admits a dual, then 
the set of admissible coding sequences is a soflc shift (Theorem 5.8). In Section 6 we 
derive formulas for the density of the absolutely continuous invariant measure and 
the measure-theoretic entropy of the one-dimensional Gauss-type maps and their 
natural extensions. We also prove that the natural extension maps are Bernoulli 
shifts. And flnally, in Section 7 we apply results of [16 to obtain explicit for- 
mulas for invariant measure for the one-dimensional maps for some regions of the 
parameter set V. 



4 



SVETLANA KATOK AND ILIE UGARCOVICI 



2. Geodesic flow on the modular surface and its representation as a 
special flow over a symbolic dynamical system 

Let 1-L = {z = X -\- iy : y > 0}be the upper half-plane endowed with the 
hyperbolic metric, T = {z e H : \z\ > 1, |Rez|<^}be the standard fundamental 
region of the modular group PSL{2,Z) = 5^(2, Z)/{±/}, and M = PSL{2,Z)\H 
be the modular surface. Let Sl-L denote the unit tangent bundle of %. We will use 
the coordinates v = (2:, Q on Sl-L^ where z G C ^ ^5 ICI = (z)- The quotient 
space PSL{2^1j)\S1-L can be identified with the unit tangent bundle of M, SM^ 
although the structure of the fibered bundle has singularities at the elliptic fixed 
points (see [11, §3.6] for details). Recall that geodesies in this model are half-circles 
or vertical half-rays. The geodesic flow {(p^} on 1-L is deflned as an R-action on the 
unit tangent bundle Sl-L which moves a tangent vector along the geodesic defined by 
this vector with unit speed. The geodesic fiow {(p^} on li descents to the geodesic 
flow {^^} on the factor M via the canonical projection 

(2.1) TT'.Sn^SM 

of the unit tangent bundles. Geodesies on M are orbits of the geodesic fiow {^p*}. 

A cross-section C for the geodesic flow is a subset of the unit tangent bundle 
SM visited by (almost) every geodesic infinitely often both in the future and in the 
past. In other words, every v e C defines an oriented geodesic ^{v) on M which 
will return to C infinitely often. The "ceiling" function g : C ^ giving the time 
of the first return to C is defined as follows: if v G C and t is the time of the first 
return ofj{v) to C, then g{v) = t. The map R:C defined by R{v) = (p^^'^^v) 
is called the first return map. Thus {(p^} can be represented as a special flow on 
the space 

= {{v, s) :veC, 0<s< g{v)}, 
given by the formula (p*{v^ s) = {v^ s-\-t) with the identification {v, g{v)) — 0)- 




Figure 2. Geodesic fiow is a special fiow 

Let A/" be a finite or countable alphabet, JV^ = {x = {ni}i^z \ ^ A/"} be the 
space of all bi-infinite sequences endowed with the Tikhonov (product) topology, 

cr : A/"^ ]\[^ defined by {(7x)i = n^+i 

be the left shift map, and A C be a closed a-invariant subset. Then (A, a) 
is called a symbolic dynamical system. There are some important classes of such 
dynamical systems. The space (JV^^a) is called the full shift (or the topological 
Bernoulli shift). If the space A is given by a set of simple transition rules which 
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can be described with the help of a matrix consisting of zeros and ones, we say 
that (A, (j) is a one-step topological Markov chain or simply a topological Markov 
chain (also called a subshift of finite type). A factor of a topological Markov chain 
is called a sofic shift;. (See [10, §1.9] for the definitions.) 

In order to represent the geodesic fiow as a special fiow over a symbolic dynamical 
system, one needs to choose an appropriate cross-section C and code it, i.e. to find 
an appropriate symbolic dynamical system (A, a) and a continuous surjective map 
Cod : A ^ C (in some cases the actual domain of Cod is A except a finite or 
countable set of excluded sequences) defined such that the diagram 

A A 

Codj |cod 

c — ^ c 

is commutative. We can then talk about coding sequences for geodesies defined 
up to a shift which corresponds to a return of the geodesic to the cross-section C. 
Notice that usually the coding map is not injective but only finite-to-one (see e.g. 
§3.2 and §5]). 

There are two essentially different methods of coding geodesies on surfaces of 
constant negative curvature. The geometric code, with respect to a given funda- 
mental region, is obtained by a construction universal for all Fuchsian groups. The 
second method is specific for the modular group and is of arithmetic nature: it uses 
continued fraction expansions of the end points of the geodesic at infinity and a 
so-called reduction theory (see [151 [14] foi" the three classical cases). Here we will 
describe a general method of arithmetic coding via (a, 6)-continued fractions that 
is based on study of the attr actor of the associated natural extension map. This 
approach, coupled with ideas of Bowen and Series [7J, may be useful for coding of 
geodesies on quotients by general Fuchsian groups. 

3. The reduction procedure 

In what follows we will denote the end points of geodesies on H by and 
and whenever we refer to such geodesies, we use {u, w) as their coordinates on 
{u ^ w). 

The reduction procedure for coding symbolically the geodesic fiow on the mod- 
ular surface via continued fraction expansions was presented for the three classical 
cases in [Mj; for a survey on symbolic dynamics of the geodesic fiow see also [15] . 
Here we describe the reduction procedure for (a, 6) -continued fractions and explain 
how it can be used for coding purposes. 

Let 7 be an arbitrary geodesic on 1-L from u io w (irrational end points), and 
w = [no, ni, . . . ]a,6. We construct the sequence of real pairs {{uk^Wk)} {k > 0) 
defined by 

(3.1) uq = wq = w and Wk-\-i = ST~^^Wk , Uk+i = ST~^^Uk . 

Each geodesic from to is P5'L ( 2, Z) -equivalent to 7 by construction. It 
is convenient to describe this procedure using the reduction map that combines the 
appropriate iterate of the map Fa^b- 

Ra,b : \ A ^ R2 \ A 
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given by the formula Ra,b{u,w) = {ST '^u^ST '^w)^ where n is the first digit in 
the (a, 6)-expansion of w. Notice that (ii^, Wk) = bi^^ 

Definition 3.1. A geodesic in 1-L from to is called (a, h) -reduced if {u^ w) G Aa,6, 
where 

Kb = Fa,b{Da,b n{a<w<b})= S{Da,b H {a < w < b}). 

According to Theorem for (almost) every geodesic 7 from u to w in 
the above algorithm produces in finitely many steps an (a, 6)-reduced geodesic 
PSL{2, Z)-equivalent to 7, and an application of this algorithm to an (a, 6)-reduced 
geodesic produces another (a, 6)-reduced geodesic. In other words, there exists a 
positive integer £ such that R^^^{u^w) G Aa,6 and Ra^b ' ^a,b ^a,b is bijective 
(with the exception of some segments of the boundary of Aa^b and their images). 

Let 7 be a reduced geodesic with the repelling point u ^ and the attracting 
point 

(3.2) W = L^0,^l,...la,6. 

Then, by successive applications of the map Ra,b^ we obtain a sequence of real pairs 



{{uk^Wk)} {k > 0) (see (3.1) above) such that each geodesic from Uk to Wk is 



(a, 6)-reduced. Using the bijectivity of the map Ra,b^ we extend the sequence (3.2) 
to the past to obtain a bi-infinite sequence of integers 

(3.3) h] = n_2, n_i, no, ni, 712, . . . 1 , 

called the coding sequence of 7, as follows. There exists an integer n_i 7^ and a real 
pair (u-i^W-i) G Aa^b such that ST~^-^W-i = w = wq and ST~^-^U-i = u = uq. 
Notice that [^^;-l]a,6 = ^-i- By uniqueness of the (a, 6)-expansion, we conclude 
that W-i = [n_i,no,ni, . . .]a,6- Continuing inductively, we define the sequence of 
integers U-k and the real pairs {u-k^W-k) G A^^^ {k > 2), where 

by ST~'^-^W-k = '^-(/c-i) and ST~^-^U-k = '^-(fe-i)- We also associate to 7 a 
bi-infinite sequence of (a, 6)-reduced geodesies 

(3.4) (...,7_2,7-i,7o,7i,72,...). 
where is the geodesic from Uk to Wk. 

Remark 3.2. Notice that all "intermediate" geodesies T~^jf^ (1 < s < rik) obtained 
from 7/c using the map Fa^b are not (a, 6)-reduced. 

Proposition 3.3. A formal minus continued fraction comprised from the digits of 
the ''past'' of (S^ , 

n_i = (n_i, n_2, n_3, . . . ) 

n_2 — 



n-3 



converges to 1/u. 

Proof By [T3l Lemma 1.1], it will be sufficient to check that \n-k\ = 1 implies 
n_/e • n_(/c+i) < 0, i.e. the digit 1 must be followed by a negative integer and the 
digit —1 must be followed by a positive integer. We use the following properties of 



APPLICATIONS OF (a, 6)-CONTINUED FRACTION TRANSFORMATIONS 



7 



the set Aa,6 that can be derived from our knowledge of the shape of the set Da,b 
determined in [16j Lemmas 5.6, 5.10, 5.11]. The upper part of Aa^b is contained in 
the region 



(3.5) 



-1,0] X 
-1,0] X 



, +00 
, +00 



U [0, 1] 



1 



, +00 



The lower part of Aa^b is contained in the region 

1 



(3.6) 



[-1,0] X 
[0, 1] X 



1 



U [0, 1] 



-oo, 



-oo, - 



if 6< 1 
if 6> 1. 

if a > -1 
if a < -1. 



Recall that (ii_(/g+i), iL'_(/e+i)) = {T'^-^^+'^^ Su-k-,T^~^^^^^ Sw-k) foi" an appropriate 
integer n_(/e+i) ^ 0. Suppose ri-k = 1. Then W-k > 0. If U-k < 0, then Su-k > 
and Sw-k < 0, and it takes a negative power of T to bring it back to (the lower 
component of) Aa^b^ i-e. n_(/c+i) < 0. The case U-k > 0, according to (3.5), can 
only take place if 6 < 1. In this case, —1/(6 — 1) < W-k < 6 + 1, which is equivalent 
to 6 > 1, a contradiction. Therefore n_/e = 1 implies n_ < 0. A similar 
argument shows that ri-k = — 1 implies n_(fc+i) > 0. We conclude that the formal 
minus continued fraction converges. In order to prove that the limit is equal to 1/u 
we use the recursive definition of the digits n_i, n_2, . . . , to write 



1 

u 



1 



n_i — U-i = n_i 



(n_i,n_ 



n_2 - U-2 

and the conclusion follows since the formal minus continued fraction converges. 
Let 

C = {zen \ \z\ = l, -1< Re z<l} 
be the upper-half of the unit circle, and 

1 



□ 



and 



C" = {z G H I |z + 1| = 1, -- < Re z < 0} 



= {zen\\z-l\ = l, < Re z < ^} 



be the images of the two vertical boundary components of the fundamental region 
T under S (see Figure [s]). 

Proposition 3.4. Every {a^b) -reduced geodesic either intersects C or both curves 
C- andC+. 



Proo f. If a, b are such that — 1 < a < and < 6 < 1, then by properties (3.5) and 
(3.6) of the set A^^b, if (u^w) G Aa,6, then — 1 < < 1 and w > — ^ or^; < — 
and hence all (a, 6)-reduced geodesies intersect C. For the case 6 > 1 we have: if 
— 1 < < 0, then either w > > b > 1 or w < — < —1, i.e. the geodesic 
intersects C; if < < 1, then (3.5) implies that k;<— ^<a<0, thus the 
corresponding geodesic intersects C if w < —1, and it intersects first and then 
if —1 < w < 0. Similarly, for the case a < — 1 we have: if < < 1, then either 



w < 



< a < —1 or w > 



1 

"6-1 



> 1, i.e. the geodesic intersects C; if —1 < u < 0, 
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then (3.6) implies that w > > b > 0^ therefore the corresponding geodesic 
intersects C if w > 1^ and it intersects first C~ and then if < < 1. □ 



Based on Proposition |3.4| we introduce the notion of the cross-section point. It 
is either the intersection of a reduced geodesic 7 with C, or, if 7 does not intersect 
C, its first intersection with C~ U 

Now we can define a map 

ip : Aa,b Sn , (p{u, w) = {z, C) 

where z G is the cross- section point on the geodesic 7 from u io and is the 
unit vector tangent to 7 at The map is clearly injective. Composed with the 



canonical projection tt introduced in (2.1) we obtain a map 

TT O : Ka^h SM. 

Let Ca,b = TT o (p{Aa^b) C SM. This set can be described as follows: (7^,6 = 
P U Qi U Q2, where P consists of the unit vectors based on the circular boundary 
of the fundamental region pointing inward such that the corresponding geodesic 
7 on the upper half-plane H is (a, 6)-reduced, Qi consists of the unit vectors based 
on the right vertical boundary of pointing inward such that either ^7 or TSj 
is (a, 6)-reduced (notice that they cannot both be reduced), and Q2 consists of the 
unit vectors based on the left vertical boundary of T pointing inward such that 
either S^y or T~^Sj is (a, 6)-reduced (see Figure [3|. Then a.e. orbit of {cp*} returns 
to i.e. Ca^b is a cross-section for {^p*}^ and Aa^b is a parametrization of Ca,b- 
The map tt o is injective, as follows from Remark 3.2 only one of the geodesies 
7, 5*7, T~^S'7, and TSj can be reduced. 





Figure 3. The cross-section (left) and its Aa^b parametrization (right) 



4. Symbolic coding of the geodesic flow via (a, 6)-continued 

FRACTIONS. 

If 7 is a geodesic on we denote by 7 the canonical projection of 7 on M. For 
a given geodesic on M that can be reduced in finitely many steps, we can always 
choose its lift 7 to H to be (a, 6)-reduced. 

The following theorem provides the basis for coding geodesies on the modular 
surface using (a, 6)-coding sequences. 
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Theorem 4.1. Let 7 be an (a, b) -reduced geodesic on H and 7 its projection to M. 
Then 

(1) each geodesic segment of j between successive returns to the cross-section 
Ca,b produces an (a^b) -reduced geodesic on H, and each reduced geodesic 
SL {2^1?) -equivalent to 7 is obtained this way; 

(2) the first return of 7 to the cross-section Ca,b corresponds to a left shift of 
the coding sequence ofj. 

Proof (1) By lifting a geodesic segment on M starting on Ca^b to we obtain a 
segment of a geodesic 7 on H that is reduced by the definition of the cross-section 
Ca^b' A coding sequence of 7 = 70 that connects uq to wq = [no, ^i, ... la, 6 5 

[70I = [. ..,n_2,n_i,no,ni,n2,...l, 

is obtained by extending the sequence of digits of wq to the past as explained in 
the previous section. 

Let us assume that wq > 0, i.e. uq > 1. The case wq < can be treated 
similarly. The geodesic ST~^^jo = 71 is reduced by Theorem Let zq and zi 
be the cross-section points on 70 and 71, respectively. Then z[ = T'^^Szi G 70; 
it is the intersection point of 70 with the circle \z — no\ = 1. We will show that 
the geodesic segment of 70, [^0,^1] projected to M is the segment between two 
successive returns to the cross-section Co, 5. Since ST~'^^{z[) = z\ is the cross- 
section point on 71, the geodesic segment [2^0,2:^] projected to M is between two 
returns to Ca.b- Recall that a geodesic in T consists of countably many oriented 
geodesic segments between consecutive crossings of the boundary of T that are 
obtained by the canonical projection of 70 to T . 

If is the intersection of 70 with C, there are two possibilities. First, when 70 
intersects T or 70 does not intersect T and ST~^^^ exists T through it circular 
boundary, and, second, when 70 does not intersect T and ST~^jo exists T through 
it left vertical boundary. In the first case the segments in are represented by the 
intersection with T of the following geodesies in H: T~^7o, T~^7o, . . . , T~"^o+^7o, 
either ST'^'^+^jo or T-^070, and either 70, or 6^-^70. 

Suppose that for some intermediate point z G 70, z G [^o^^i] the unit vector 
tangent to 70 at (z, (") is projected to Ca^b- By tracing the geodesic 70 inside J^, 
we see that (z, Q must be projected to (z, () with z on the boundary of and ( 
directed inward. Then the geodesic through (z, () 

(a) enters T through its vertical boundary and exits it also through the vertical 
boundary, 

(b) enters through its vertical boundary and exits through its circular bound- 
ary, or 

(c) enters T through its circular boundary and exits through its vertical bound- 
ary. 

The following assertions are implied by the analysis of the attractor Da^b- In case 
(a), T~^ST~^jo is not reduced for 1 < 5 < no since 5 < no, T~^wo > 6, hence 
ST-'wo > -I, i.e. {ST-'uo,ST-'wo) ^ Da,b^ therefore 

{T-^ST-'uo,T-^ST-'wo) ^ Aa,b- 

In case (b), either the segment T~'^^jo exits through the circular boundary of J^, 
ST~'^^jo = 7i is reduced, and we reached the point zi on the cross-section. If 
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the segment T~^°+^7o intersects the circular boundary of T ^ ST~'^^~^^jo is not 
reduced. In case (c), ST~'^^~^^ is not reduced. 

In the second case the first digit of wq^ no = 2. This is because no = 1 would 
imply b -\- 1 < w < which is impossible. Thus ST~'^jo = 71 is reduced. In 

this case the geodesic in T consists of the intersection with of a single geodesic 
S'T~^7o that enters T through its right vertical and leave it through its left vertical 
boundary, since {TS)T{ST~^jq) = ST~'^jo = 71 is reduced. In all cases the 
geodesic segment [2^0,^1] projected to M is between two consecutive returns to 

Ca,b' 



If zq ^ C, by Proposition |3.4[ since wq > 0, G C~ . Notice that this implies 
that a < —1 and no = 1, and 71 = 6'T~^7o is reduced. In this case the geodesic 
in also consists of the intersection with of a single geodesic Sjo that enters 

through its right vertical and leave it through its left vertical boundary, since 
{TS)T{Sjo) = ST~^jo = 7i is reduced, and hence the geodesic segment [^0,^1] 
projected to M is between two consecutive returns to Ca^b- Continuing this argu- 
ment by induction in both positive and negative direction, we obtain a bi-infinite 
sequence of points 

(. . . , Z_2, 2^0, 2^1, ^2, . . . ), 

where is the cross-section point of the reduced geodesic in the sequence of 70, 
that represents the sequence of all successive returns of the geodesic 70 in M to the 
cross-section Ca,b- 

If 7o is a reduced geodesic in 5'L ( 2, Z) -equivalent to 70, then both project to 
the same geodesic on M. Therefore, the cross-section point zq of 70 projects on 
Ca^b to a cross-section point z^ of 7^ for some k. This completes the proof of (1). 

(2) Since 71 = 6^-^070, wi = ST'^'^wo = L^i, n2, . . . la,6. The first digit of 
the past is evidently no, and the remaining digits are the same as for 70. Thus (2) 
follows. □ 

The following corollary is immediate. 

Corollary 4.2. Ifj^ is SL (2^1?) -equivalent to ^, and both geodesies can be reduced 
in finitely many steps, then the coding sequences of 7 and 7' differ by a shift. 

It implies a very important property of (a, 6)-continued fractions that escapes a 
direct proof. 

Corollary 4.3. (The Tail Property) For almost every pair of real numbers 
that are SL (2^ I?) -equivalent, the ''tails" of their (a, b)- continued fraction expansions 
coincide. 



Remark 4.4. The set of exceptions in Corollary |4.3| is the same as the one described 
m Theorem [TlJS). 

Thus we can talk about coding sequences of geodesies on M. To any geodesic 7 



that can be reduced in finitely many steps we associate the coding sequence (3.3) 



of a reduced geodesic 6'L ( 2, Z) -equivalent to it. Corollary 4.2 implies that this 
definition does not depend on the choice of a particular representative: sequences 
for equivalent reduced geodesies differ by a shift. 

Let Xa^b be the closure of the set of admissible sequences and a be the left shift 
map. The coding map Cod : Xa^b Ca,b is defined by 

(4.1) Cod([. . . , n_2, n_i, no, ni, ...]) = (1/ (n_i, n_2, [uq, ni, . . .]a,b)' 
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This map is essentially bijective. 

The symbolic system (Xa,6,cr) C (A/'^,cr) is defined on the infinite alphabet 
A/" C Z \ {0}. The product topology on A/"^ is induced by the distance function 

m 

where x = {ni)^x' = {n[) G AT^, and m = max{/c | = for \i\ < k}. 
Proposition 4.5. The map Cod is continuous. 

Proof. If x^) < ^, then the (a, 6)-expansions of the attracting end points w{x) 
and w{x') of the corresponding geodesies given by ( |3.2[ ) have the same first m 
digits. Hence the first m convergents of their (a, 6)-expansions are the same, and 
using the properties of (a, b) continued fraction and the rate of convergence of[16l 
Theorem 2.1] we obtain \w{x) — w{x')\ < ^. Similarly, the first m digits in the 
convergent formal minus continued fraction of and —r-r\ are the same, and 

u{x) u{x') ' 

hence \u{x) — u{x')\ < (^)l < 1.^ Therefore the geodesies are uniformly ~ 

close. But the tangent vectors v{x)^ v{x^) G Ca,b are determined by the intersection 
of the corresponding geodesic with the unit circle or the curves and C~ . Hence, 
by making m large enough we can make v{x') as close to v{x) as we wish. □ 

In conclusion, the geodesic fiow becomes a special fiow over a symbolic dynamical 
system (Xa,6, a) on the infinite alphabet A/" C Z \ {0}. The ceiling function ga,b{x) 
on Xa^b coincides with the time of the first return of the associated geodesic 7(0;) 
to the cross-section Ca^b- One can establish an explicit formula for ga,b{x) as the 
function of the end points of the corresponding geodesic 7(x), u{x)^ w{x)^ following 
the ideas explained in [8 . If — 1 < a < and < 6 < 1, then ga,b{x) is cohomologous 
to 2 log more precisely, 

ga,b{x) = 2\og\w{x)\^\ogh{x)-\ogh{(jx) where h{x) = ' ^ \ ^ /'^ . 

1 — u[x)'^ 

5. Dual codes 



We have seen that a coding sequence for a reduced geodesic from to (3.3) is 
comprised from the sequence of digits in (a, 6)-expansion of w and the "past" , an 
infinite sequence of non-zero integers, each digit of which depends on w and u. In 
some special cases the "past" only depends on and, in fact, it will coincide with 
the sequence of digits of by using a so-called dual expansion to (a, b). 

Let y) = (— —x) be the refiection of the plane about the line y = —x. 

Definition 5.1. If ip{Da,b) coincides with the attractor set Da'^b' for some (a', b') G 
7^, then the (a^, 6')-expansion is called the c/i^a/ expansion to (a, b). If {a\ b') = (a, 6), 
then the (a, 6)-expansion is called self-dual 

Example 5.2. The classical situations of ( — 1,0)- and ( — 1, l)-expansions are self- 
dual. Two more sophisticated examples ( ^~2^ ) §)' respectively, 
are shown in Figure |4j 

Example 5.3. The expansions ( — ^,1 — ti > 1, satisfy a weak cycle property 
and have dual expansions that are periodic. A classical example in this series is the 
Hurwitz case (— ^, ^) whose dual is (^— ~^~^^ ) (see [9l[T3]). Their domains are 
shown in Figure P 
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1 





Figure 4. Domains of self-dual expansions 




(a,6) = (-i,i) 



(a',6') = (^,-4^) 



Figure 5. Dual expansions 



The following result gives equivalent characterizations for an expansion to admit 
a dual. 

Proposition 5.4. The following are equivalent: 

(i) the {a, b) -expansion has a dual; 

(ii) the boundary of the lower part of the set Da^b does not have y-levels with 
a < y < 0, and the boundary of the upper part of the set Da,b does not have 
y-levels with < y < b; 

(iii) a and b do not have the strong cycle property. 

Proof If the (a, 6)-expansion has a dual (a', 6')-expansion, then the parameters a', b' 
are obtained from the boundary of Da^b as follows: the right vertical boundary of 
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the upper part of Da,b is the ray x = 1 — b\ and the left vertical boundary of the 
lower part of Da^b is the ray x = — 1 — a^ Now assume that (ii) does not hold. 
Then at least one of the parameters a, b has the strong cycle property, and either 
the left boundary of the upper part of Aa,6 or the right boundary of the lower part 
of Att,6 is not a straight line. Assume the former. Then the reflection of Da^b with 
respect to the line y = —x is not Da'^b' since the map Fa'^b' is not bijective on it: 
the black rectangle in Figure [g] belongs to it, but its image under T~^, colored in 
grey, does not. Thus (i)^(ii). 




Figure 6. Dual expansions and Da^b 

Conversely, let the vertical line x = 1 — b' he the right boundary of the upper 
part of Da^b and the vertical line x = —1 — he the left boundary of the lower 
part of Da,b' Let [x^, oo] x {a} he the intersection of Da,b with the horizontal line 
at the level a, and [— 00,^5] x {b} he the intersection of Da^b with the horizontal 
line at the level b. Then = ^ and b' = We also see that 1 — b' = —7, 
where t = Xb or t < Xb if [t^Xb] x {0} is a segment of the boundary of Da^b- Then 
—b^ + 1 = — ^ < which implies b' — a' > 1. By Lemma 5.6 of [TF Xb < —1 and 

> 1, therefore 

(5.1) - 1 < a' < < 6' < 1, 
and 

(5.2) Aa,b = Da^b n {{u, w) eR'^ : -b' <u< -a'}. 

We now show that ip^Da^b) = Da'^b' is the attractor for Fa' ^b' ^ where 

(5.3) Fa'^b'=^oF-lo^-\ 

For (u^w) G Da'^b' with a' < w < b' we have iIj~^{u^w) = {—w^—u) with —b' < 



u < —a'^ so 'il)~^(u^w) G Aa,6 by (5.2), hence F~l{—w^—u) = {l/w^l/u)^ and 
Fa'^b'{u^w) — {—1/u^—l/w). For {u^w) G Da'^b' with w > b' we have il)~^{u^w) = 
(— k;, —u) with u < —b\ so F~^{—w, —u) = (— il' + I, — ii+1), and Fa' ^b' {u^ w) = {u — 
l^w — 1). Similarly, for (u^w) G Da'^b' with < a' we have tlj~^{u^w) = {—w^ —u) 
with u > —a'^ so F~l{—w^ —u) = — —u — 1)^ and Fa'^b'i^^ w) = {u-\-l^w-\-l). 
This proves that (ii)=>(i). 
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Notice that (ii) and (iii) are equivalent by Theorems 4.2 and 4.5 of [16]. □ 

Remark 5.5. Notice that if an (a, 6)-expansion has a dual, then — l<a<0<6<l. 
This follows from ( |5.1[ ) and the fact that the relation of duality is symmetric. 

Theorem 5.6. If an {a ^b)- expansion admits a dual expansion {a'^h'), and is an 
(a^b) -reduced geodesic, then its coding sequence 

(5.4) [70] = [. ..,n_2,n_i,no,ni,n2,...l, 

is obtained by juxtaposing the {a^b) -expansion of wq = [no, ni, n2, . . . ]a,6 ci'^d the 
{a' ^b')- expansion ofl/uo = [n_i, n_2, . . . ]a',6' • This property is preserved under 
the left shift of the sequence. 

Proof We will show that the digits in the (a^, 6')-expansion of I/uq coincide with 
the digits of the "past" of (5.4). By (5.3), the following diagram 



Aa,6 - 


Sip 


Aa 


',b' 








Ra 


Kb - 






',b' 



is commutative. The pair {uq^wq) G A^^^, therefore {Suq^Swq) G SAa^b C Da^b^ 
and {l/wo,l/uo) G Aa'^b'- The first digit of the (a^, 6')-expansion of 1/iio is so 

(1/^0,1/^0) = {ST-^-^{l/wo)^ST-^-^{l/uo)) 

maps Aa'^b' to itself. Then 

:= R-l{uo,wo) = (T"-^5iXo,T"-^^^o) e A,,5 

and {ST~'^-^U-i^ ST~'^-^W-i) = {uq^wq). Also W-i = [n_i, no, ni, . . . ]a,6, and 

ST-^-^{l/uo) = 1/ii-l = L^-2,...la',6- 

Continuing by induction, one proves that all digits of the "past" of the sequence 
(5.4) are the digits of the (a^, 6')-expansion of I/uq. 

In order to see what happens under a left shift, we reverse the diagram to obtain: 



Aa,6 - 






,b' 






1 




Aa,6 - 






,b' 



Since the first digit of (a, 6)-expansion of wq is no, 

maps Aa,6 to itself. Then (ui^wi) := {ST~^^uo^ ST~^^wo) and wi = [ni, n2, . . • la,6- 
Also 

{l/wr,l/m)=R-\^{l/woA/uo) = (To^^(l/^o),To^^(l/^io)), 
hence 1/ui = [no, n_i, n_2, . . . □ 



Remark 5.7. Under conditions of Theorem |5.6[ if 70 projects to a closed geodesic 
on Tkf, then its coding sequence is periodic, and wq = [no,ni, . . . , n^ 
[um, ' . . ,ni,no] a',b'- 

Theorem 5.8. // an {a, b)- expansion admits a dual expansion, then the symbolic 
space (Xa,6,cr) is a sofic shift. 
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Proof. The "natural" (topological) partition of the set Aa,6 related to the alphabet 
JV is Aa^b = UnGA^An, where are labeled by the symbols of the alphabet JV and 
are defined by the following condition: = {{u^w) G A^^^ | no{u,w) = no{w) = 
n}. In order to prove that the space (X^^b, a) is sofic one needs to find a topological 
Markov chain {Ma.b^ ^) and a surjective continuous map h : Ma^h ^ ^a,b such that 
h o r = a o h. 

Notice that the elements A^ are rectangles for large n; in fact, at most two 
elements in the upper part and at most two elements in the lower part of A^^^ are 
incomplete rectangles (see Figure [7|. 





Figure 7. The partition of Aa^b and its image through Ra^b- 



Since Aa^b has finite rectangular structure, we can sub-divide horizontally these 
incomplete rectangles into rectangles, and extend the alphabet JV by adding sub- 
scripts to the corresponding elements of JV. For example, if A2 is subdivided into 
two rectangles, A2 = U^^iA2., the "digit" 2 will give rise to two digits, 2i,22 in 
the extended alphabet A/"' (see Figure [t]). We denote the new partition of A^^^ by 
UnGAT'^n- Noticc that it consists of rectangles with horizontal and vertical sides. 
Since the first return R to Aa^b corresponds to the left shift of the coding sequence 
X associated to the geodesic (u^w)^ we see that x = {^fej^oc where Uk is defined 
by R^{u,w) G A^^. Now we define the symbolic space Ma^b as follows: to each 
sequence x G Xa^b we associate a geodesic (u^w) by (4.1), and define a new coding 
sequence y = {^/cj^oc where rrik is defined by R^{u^w) G M^^, and r is the left 
shift. 
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We will prove that {Ma^b^'^) is a topological Markov chain. For this, in accor- 
dance to [2 , Theorem 7.9], it is sufficient to prove that for any pair of distinct symbol 
n, m G JV' ^ R{Mn) and Mm either do not intersect, or intersect "transversally" i.e. 
their intersection is a rectangle with two horizontal sides belonging to the horizon- 
tal boundary of Mm and two vertical sides belonging to the vertical boundary of 
R{Mn). Let us recall that — l<a<0<6<l (see Remark 5.5). Therefore, 
if Mn = is a complete rectangle, it is, in fact, a 1 x 1 square, and its image 
under R is an infinite vertical rectangle intersecting all Mm transversally. If Mn is 
obtained by subdivision of some and belongs to the lower part of A^^^, its hori- 
zontal boundaries are the levels of the step-function defining the lower component 
of Da,b^ and by Proposition 5.4, since the lower boundary of Da^b does not have 
^/-levels with a < < 0, its image is a vertical rectangle intersecting only the lower 
component of Da^b whose horizontal boundaries are the levels of the step-function 
defining the lower component of Da^b- Therefore, all possible intersections with Mm 
are transversal. A similar argument applies to the case when belongs to the 
upper part of Aa^b- The map h : Ma^b ^a,b is obviously continuous, surjective, 
and, in addition, h o r = a o h. □ 



6. Invariant measures and ergodic properties 

Based on the finite rectangular geometric structure of the domain Da^b and the 
connections with the geodesic flow on the modular surface, we study some of the 
measure-theoretic properties of the Gauss- type map fa^b • ^) 



(6.1) f\b{x) 



1 


1" 




X 


X 


a,b 



/a,6(0) =0. 



Notice that the associated natural extension map Fa^b 

(6.2) Fa^b{x, y) = ( /a,6(^), \ \ / 1 

is obtained from the map Fa^b induced on the set Aa^b by the change of coordinates 

(6.3) x = —l/w^ y = u 

(or, equivalently, on the set Da^b H {{u^ w)\a < w < b} by the change of coordinates 
X = w^y = —1/u). Therefore the domain Aa^b of Fa^b is easily identified knowing 
Aa^b and may be considered as its "compactification" . 

Many of the measure-theoretic properties of fa^b and Fa,6 (existence of an ab- 
solutely continuous invariant measure, ergodicity) follow from the fact that the 
geodesic flow cp^ on the modular surface M can be represented as a special flow 
{Ra,b,K,b,ga,b) on the space 

K^,b = {(^^^'^) • (^'^) ^ ^a,6, < t < ga^b{u,w)} 

(see Section 2). We recall that Ra^b = Fa,b\Aa,b 9a,b is the ceiling function (the 

time of the first return to the cross-section Ca^b) parametrized by {u,w) G A^^^. 

We start with the fact that the geodesic flow {cp^} preserves the smooth (Liou- 

11 N 7 dudwdt , x 1 1111 

ville) measure dm — ttt (see, e.g., 3 ), hence Rab preserves the absolutely 

[w — uy 
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_ dudw . ^ ^ _ _ r— ^ 

continuous measure dp = -tz. Using the change oi coordinates (6.3), the map 

[w — u)"^ ^ ' 

A 1111. 7 dxdy 

hah preserves the absolutely continuous measure dv = -77. 

(1 + x?/)^ 

The set Ka^h has finite measure dp if a / and 6^0, since it is uniformly 



bounded away from the line A = {(ii,^;) \ u = w} dW' (see relations (3.5) and 



(3.6)). In this situation, we can normalize the measure dp to obtain the smooth 



probability measure 

dp dudw 



(6.4) dpa^l 



Ka,h Ka^b{w-uy 



where Ka^b = p{^a,b)- Similarly, if a 7^ and 6 7^ 0, the map Fa^b preserves the 
smooth probability measure 

(6.5) dUa ^"^^ 



KaA^^xyy 



and Ka,b = p{K,b) = ^{K.h)- 

Returning to the Gauss- type map, fa^bi one can obtain explicitly a Lebesgue 
equivalent invariant probability measure pa,h by projecting the measure Va,h onto 
the x-coordinate (push- forward); this is equivalent to integrating Va,h over Ka,h with 
respect to the ^/-coordinate as explained in [4]. 

We can immediately conclude that the systems (i^a, 65^0,6) and {fa,b^Pa,b) are 
ergodic from the fact that the geodesic fiow {cp^} is ergodic with respect to dm. 
By using some well-known results about one dimensional maps that are piecewise 
monotone and expanding, and the implications for their natural extension maps, 
we can establish stronger measure-theoretic properties: (/a, 6, Pa, 5) is exact, and 
{Fa^bi^a,b) is a Bernoulli shift. Here we follow the presentation from [23] based on 

Eu'lTH].' 

Theorem 6.1. For any a 7^ and b ^ 0, the system {fa,b^Pa,b) 'is exact and its 
natural extension {Fa^b^^a,b) 'is a Bernoulli shift. 

Proof. Let us consider first the case — l<a<0<6<l. The interval (a, b) admits 
a countable partition ^ = {-^i}iGZ\{o} of open intervals and the map /a, 5 satisfies 
conditions (A), (F), (U) listed in [25]. Condition (A) is Adler's distortion estimate: 

(^) ■■ f'lJUib? is bounded on X = Uiez\{o}^«, 
condition (F) requires the finite image property of the partition ^, 

{F) : fa,b{0 = {fa,b{^i)}iez\{o} is finite, 
while condition (U) is a uniformly expanding condition 

(U): |/;,|>T>lonX 

Let m > and n > be such that a — m < —l/b<a — m — l and 6 + n < —1/a < 
6 + n + 1. Consider the open intervals 



^1 = 7 ,b\, Xi={ : , — - for i > 2 

V a — m — 1 I \ a — m — i a — m — i + 1 
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and 

\ o + n + ly \ o + n + z — 1 b-\-n-\-iJ 

The map fa^b satisfies conditions (A), (F), (U) with respect to the partition ^ = 
{Xi}iez\{o}' Indeed, \f'lb/Ua,hf\ < 2 on X, the cohection of images /a,6(0 consists 
of four sets /a,6(^i), /a,6(^-i), " 1,^), («, « + 1), and |/^ J > min{^, ^} > 1 
on X. Zweimiiller ^] showed that any one-dimensional map for which conditions 
(A), (F), (U) hold is exact and satisfies Rychlik's conditions described in [18 , hence 
its natural extension map is Bernoulli. 

We analyze now the case h > 1. Let K > {) he the smallest integer such that 
6(a + 1)^ < 1. We will show that there exists 7 > 1 such that, for every x G 
riz^o /a"6(^)' some iterate f^^^{x) with n < + 1 is expanding, i.e. \{fa^^)\x)\ > 7. 
(For the rest of the proof, we simplify the notations and let / denote the map /a, 5.) 
Notice that if x G nr=o /~*(^)' then is differentiable at x and 



dx {xf{x)--- f^-^{x)Y 

Assume that ah > —1. We look at the following cases: 

(i) If a < X < 0, then b-l< f{x) < 6, and \xf{x)\ < \ab\ < 1. 

(ii) If < X < 6, then a < f{x) < a + 1. Let iiT be such that b{a + 1)^ < 1. 
Then either there exists 1 < n < K such that < P{x) < a + 1 for z = 
1, 2, . . . , n - 1 and a < /^(x) < 0, or < /'{x) < a + 1 for z = 1, 2, . . . , i^. 
In the former case we have that 

(6.6) \xf{x)---rix)\ < \ab{a+ir-'\ < 1, 
while in the latter case 

(6.7) |x/(x)---/^(x)|<|6(a + l)^|<l. 

In the case a6 = — 1, let r, e > be sufficiently small such that 

b < -l/(a + r) < 6 + 1 and a - 1 < -l/{b - e) < a . 

We have: 

(i) If a < X < a-\-r^ then b—l< f{x) < — l/(a+r), and \xf{x)\ < \a/ (a+r)| < 
1. If a + r < X < 0, then \xf{x)\ < \b{a ^ r)\ < 1. 



(ii) If b — e < X < b, then < f{x) < a + 1 and one has either (6.6) with n > 2 



or (6.7). If < X < 6 — e, then one has (6.6) or (6.7) where b is replaced by 
b-e. 

In conclusion, there exists a constant 7 > 1 such that for every x G Hilo faH^) 
some iterate f^bi^) with n < K -\- 1 satisfies the condition K/^ 6)^(^)1 ^ 7- This 
implies that the iterate /j^^, with N — {K + 1)!, is uniformly expanding, i.e. it 
satisfies property (U). Since properties (A) and (F) are automatically satisfied by 
any iterate of fa;b (see [23]), we have that F^^ is Bernoulli. Using one of Ornstein's 



results [T71 Theorem 4, p. 39], it follows that Fa^h is Bernoulli. □ 
The next result gives a formula of the measure theoretic entropy of (Fa, 6 7 ^0,5)- 
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Theorem 6.2. The measure-theoretic entropy of {Fa^b^^a,b) '^s given by 

1 TT^ 



(6.8) h,.AFa,b)=^ . 

J^a,b o 

Proof. To compute the entropy of this two-dimensional map, we use Abramov's 
formula 



where rh is the normalized Liouville measure dm = — ^ttv. It is well-known that 
m{SM) = 7r^/3 (see [3]) and /im({0*}) = 1 (see, e.g., The measure dm can be 

represented by the Ambrose-Kakutani theorem [5^ as a smooth probability measure 
on the space A^''^^ 

(6.9) dm = 



/a ^9a,bdpa,l 



where dpa^b is the probability measure on the cross-section Aa^b given by (6.4). This 
implies that 

_ ^ dpdt dm 
dm = = 

Ka,b /a ^9a,bdpa,b Ka^b /a ^ 9a,bdpa,b 

Therefore Ka^b Xv^ , Qa.bdpa.b = m{SM) = tt'^/S and 

1 TT^ 



b{Fa,b) = hp^ ^^{Ra,b) = / ga,bdpa,b j. ^ 
J^a,b ^ci,b 3 



□ 



Since {Fa^b^i^a^b) is the natural extension of {fa,b^Pa,b)^ the measure-theoretic 
entropies of the two systems coincide, hence 

(6.10) h^.Akb) = j^"^ ■ 

As an immediate consequence of the above entropy formula we derive a growth 
rate relation for the denominators of the partial quotients Pn/Qn of {cl-, &)-continued 
fraction expansions, similar to the classical cases. 

Proposition 6.3. Let {qn{x)} he the sequence of the denominators of the partial 
quotients Pn/ qn associated to the {a ^h)- continued fraction expansion of x e [a, 6). 
Then 

(6.11) lim = -hp^ ,{fa,b) = T^-^ for a.e. x. 

n^oo n 2 J^a,b ^ 

Proof The proof is similar to the classical case: using the Birkhoff 's ergodic theo- 
rem one has 

b 



\ogqn{x) f 
hm = - / log \x\dpa,b • 



At the same time, Rokhlin's formula tells us that 

pb nb 

hp^^,{fa,b)= log|/^5|d/ia,5 = -2 / \og\x\dpa,b^ 
J a J a 

hence the conclusion. □ 
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7. Some explicit formulas for the invariant measure /ia,6 



In order to obtain explicit formulas for jUa^b and h^^ i,{fa,b)^ one obviously needs 
an explicit description of the domain Da^b- In [HI we describe an algorithmic 
approach for finding the boundaries of Da^b for all parameter pairs (a, b) outside of 
a negligible exceptional parameter set E. Let us point out that the set Da^b niay 
have an arbitrary large number of horizontal boundary segments. The qualitative 
structure of Da^b is given by the cycle properties of a and b. This structure remains 
unchanged for all pairs (a, b) having cycles with similar combinatorial complexity. 
For a large part of the parameter set the cycle descriptions are relatively simple 
(see Section 4]) and we discuss it herein. 

In what follows, we focus our attention on the situation — l<a<0<6<l, 
and due to the symmetry of the parameter set with respect to the parameter line 
a = —b we assume that a < —b. 

We treat the case I < < b -\- 1 and a < — < a + 1 (for some m > 
1). The coordinates of the corners of the boundary segments in the upper region 
Da,b n {{u, w)\u < 0, a < w < b} are given by 

(-2,6-1), (-|t-25(6-1)),...,(-^,(T-25)(-i)(6-1)), (-1,-1 -l) 

while the corners of the boundary segments in the lower region Da^b H {{u^w)\u > 
O^a < w <b} are given by 



(^m, + m^, (m + l,a + 1) . 



Therefore the set Aa^b is given by 

m—l 



Kb = U [iT-'sr-\b - 1), {T-^srib - 1)] X [0, 

U[{T-^Sr-\b-l),---l] X [0,^^]u[---l,6] X [0,1] 

CL TTXi I X CL 

U [a,-^ +m] X [- — ,0] U [-^ +m,a + l] X [ ^,0] 

b m b m + 1 



Theorem 7.1. //l<-^<6 + l and a<-\^m<a^l, then 

1 



-ha^b{x)dx, 



where Ka^b = log[(m — a)(l + 6)^ ^] and ha^b{x) = h'^bi^) + bi^) ^^^^ 
1 



Kbi"^) 



1 



m+l 



if (T-'^Sy-^ib -l)<x< (T-'^Syib - = 1, . . . , m - 1 
if {T-^S)'^-\b -l)<x<-^-l 
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Figure 8. Typical domain Aa,6 for the case studied 



and 



1 



m — X 
1 



ifa<x<—j--\-m 



if 



m < X < a -\- 1 . 



- m + 1 — X ^ 
Proof. The density formulas are obtained from the simple integration result 

(7.2) ^^1.1/1 I \ d 



(1 + xy)'^ X \1 -\- dx 1 -\- cx J 1 -\- dx 1 -\- cx 



For the density in the upper part of Aa,6, y > 0, all integrals have the lower 



boundary c = 0, hence the result of (7.2) becomes l/{x + 1/d). This gives the 



description of h'^^{x). For the density in the lower part of Aa^b^ V all integrals 
have the upper boundary (i = 0, hence the result —l/{ — l/c — x) and the description 
of h~^{x). By a somewhat tedious computation, we get 



/ ha,b{x)dx = log[(m - a)(l + 6)^ 



and this completes the proof. 



□ 
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